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There is growing evidence from both experiment and numerical studies that low half-odd integer 
quantum spins on a kagome lattice with predominant antiferromagnetic near neighbor interactions 
do not order magnetically or break lattice symmetries even at temperatures much lower than the 
exchange interaction strength. Moreover, there appear to be a plethora of low energy excitations, 
predominantly singlets but also spin carrying, which suggest that the putative underlying quantum 
spin liquid is a gapless "critical spin liquid" rather than a gapped spin liquid with topological order. 
Here, we develop an effective field theory approach for the spin-1/2 Heisenberg model with easy- 
plane anisotropy on the kagome lattice. By employing a vortex duality transformation, followed 
by a fermionization and flux-smearing, we obtain access to a gapless yet stable critical spin liquid 
phase, which is described by (2+l)-dimensional quantum electrodynamics (QED3) with an emergent 
SU(8) flavor symmetry. The specific heat, thermal conductivity, and dynamical structure factor are 
extracted from the effective field theory, and contrasted with other theoretical approaches to the 
kagome antiferromagnet. 



I. INTRODUCTION 

Kagome antiferromagnets are among the most extreme 
examples of frustrated spin systems realized with nearest- 
neighbor interactions. Both the frustration in each trian- 
gular unit and the rather loose "corner-sharing" aggre- 
gation of these units into the kagome lattice suppress the 
tendency to magnetically order. On the classical level, 
kagome spin systems are known to exhibit rather spe- 
cial properties: nearest-neighbor Ising and XY models 
remain disordered even in the zero-temperature limit, 
while the 0(3) system undergoes order- by-disorder into 
a coplanar spin structured 

Quantum kagome antiferromagnets, which are much 
less understood, provide a fascinating arena for the pos- 
sible realization of spin liquids. Early exact diagonal- 
ization and series expansion studies^ 3 -!^ as well as fur- 
ther exhaustive numerical works jSiLSai provide strong ev- 
idence for the absence of any magnetic order or other 
symmetry breaking in the nearest neighbor spin-1/2 sys- 
tem. Moreover, a plethora of low energy singlet excita- 
tions is found, below a small (if non-zero) spin gap. But 
the precise nature of the putative spin liquid phase in 
this model has remained elusive. 

On the experimental front, several 
mcnsional materials with magnetic 
a kagome arrangement have been 
eluding SrCi'8- a; Ga4+ a: Oig with Cr 3+ 
momenta±&ii±^^, jarosites KM 3 (OH) 6 (S0 4 )2 with 
M=Cr 3 + or Fc 3+ (S = 5/2) moments^^^, and 
volborthite Cu 3 V 2 7 (OH) 2 - 2H 2 with Cu 2 + (S = 1/2) 
moments^. Recently, herbertsmithite ZnCu3(0H)eCl2 
which also has Cu 2+ moments on a kagome lattice was 
synthesized for which no magnetic order is observed 
down to 50 mK despite the estimated exchange constant 
of 300 K-I^i 21 The second layer of 3 He absorbed 
on graphoile is also believed to realize the kagome 
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The suppression of long-range spin cor- 



relations or other signs of symmetry breaking down 
to temperatures much lower than the characteristic 
exchange energy scale is manifest in all of these ma- 
terials, consistent with expectations. Moreover, there 
is evidence for low energy excitations, both spin car- 
rying and singlets. But the ultimate zero-temperature 
spin liquid state is often masked in these systems by 
magnetic ordering or glassy behavior at the lowest 
temperatures, perhaps due to additional interactions or 
impurities, rendering the experimental study of the spin 
liquid properties problematic. New materials and other 
experimental developments are changing this situation, 
and the question of the quantum spin liquid ground 
state of the kagome antiferromagnet is becoming more 
prominent. 

There are two broad classes of spin liquids which have 
been explored theoretically, both in general terms and 
for the kagome antiferromagnet in particular. The first 
class comprise the "topological" spin liquids, which have 
a gap to all excitations and have particle- like excitations 
with fractional quantum numbers above the gap. Ar- 
guably, the simplest topological liquids are the so-called 
Z2 spin liquids, which support a vortex like excitation - 
a vison - in addition to the spin one-half spinon. For the 
kagome antiferromagnet, Sachde\*21 and more recently 
Wang and Vishwanath 2 ^ have employed a Schwinger bo- 
son approach to systematically access several different 
Zi spin liquids. For a kagome antiferromagnet with easy 
axis anisotropy and further neighbor interactions, Ba- 
lents et al2^ unambiguously established the presence of a 
Z 2 spin liquid, obtaining an exact ground state wavefunc- 
tion in a particular limit. Quantum dimer models on the 
kagome lattice can also support a Z 2 topological phased. 
Spin liquids with topological order and time reversal sym- 
metry breaking, the chiral spin liquids which are closely 
analogous to fractional quantum Hall states, have been 
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found on the kagome lattic e) 28 ' 29 within a fermionic repre- 
sentation of the spins. But all of these topological liquids 
are gapped, and cannot account for the presence of many 
low excitations found in the exact diagonalization studies 
and suggested by the experiments. 

A second class of spin liquids - the "critical" or "al- 
gebraic spin liquids" (ASL) - have gapless singlet and 
spin carrying excitations. Although these spin liquids 
share many properties with quantum critical points such 
as correlation functions falling off as power laws with non- 
trivial exponents, they are believed to be stable quantum 
phases of matter. Within a fermionic representation of 
the spins, Hastings^ - explored an algebraic spin liquid 
on the kagome lattice, and more recently Ran et al£^ 
have extended his analysis in an attempt to explain the 
observed properties of herbertsmithite ZnCu3(OH) 6 Cl2. 

In this paper, motivated by the experiments on the 
kagome materials and the numerical studies, we pursue 
yet a different approach to the possible spin liquid state 
of the kagome antiferromagnet. As detailed below, for 
a kagome antiferromagnet with easy-plane anisotropy we 
find evidence for a new critical spin liquid, an "algebraic 
vortex liquid" (AVL) phase. Earlier we had introduced 
and explored the AVL in the context of the triangular 
XXZ antiferromagnet in Refs. USUI Compared to the 
slave particle techniques for studying frustrated quantum 
antifcrromagncts, the AVL approach is less microscopi- 
cally faithful to a specific spin Hamiltonian, but has the 
virtue of being unbiased. 

Our approach requires the presence of an easy-plane 
anisotropy, which for a spin one-half system can come 
from the Dzyaloshinskii-Moriya interaction, although 
this is usually quite small. But for higher half-integer 
spin systems, spin-orbit coupling allows for a single ion 
anisotropy D(S Z ) 2 which can be appreciable. With 
D > this leads to an easy-plane spin character, and 
from a symmetry stand point our analysis should be rel- 
evant. Moreover, a very interesting unpublished exact 
diagonalization study by Sindzingre^ suggests that an 
unusual spin liquid phase is also realized for the nearest- 
neighbor quantum XY antiferromagnet. He finds a small 
gap to S z = 1 excitations, but below this gap there 
is a plethora of S z = states, which is reminiscent of 
the many singlet excitations below the triplet gap in the 
Heisenberg SU(2) spin model£ Thus, based on the exact 
diagonalization studies, the easy plane anisotropy does 
not appear to gap out the putative critical spin liquid of 
the Heisenberg model, although it will presumably mod- 
ify its detailed character. 

Our approach to easy-plane frustrated quantum anti- 
ferromagnets focuses on vortex defects in the spin con- 
figurations, rather than the spins themselves. In this 
picture, the magnetically ordered phases do not have 
vortices present in the ground state. But there will be 
gapped vortices, and these can lead to important effects 
on the spectrum. For example, vortex-antivortex (i.e. ro- 
ton) excitations can lead to minima in the structure fac- 
tor at particular wavevectors in the Brillouin zone ; 36 ' 37 ' 38 



in analogy to superfluid Hc-4. As we demonstrated for 
the triangular antiferromagnet and explore below for the 
kagome lattice, our approach predicts specific locations 
in the Brillouin zone for the roton minima, which could 
perhaps be checked by high order spin wave expansion 
techniques. On the other hand, the presence of mobile 
vortex defects in the ground state can destroy magnetic 
order, giving access to various quantum paramagnets. If 
the vortices themselves condense, the usual result is the 
breaking of lattice symmetries, such as in a spin Peierls 
state. But if the vortices remain gapless one can access 
a critical spin liquid phase. 

For frustrated spin models the usual duality transfor- 
mation to vortex degrees of freedom does not resolve the 
geometric frustration since the vortices are at finite den- 
sity. However, by binding 27r-flux to each vortex con- 
verting them into fermions coupled to a Chern-Simons 
gauge field followed by a simple flux-smearing mean-field 
treatment gives a simple way to describe the vortices. 
The long-range interaction of vortices actually works to 
our advantage here since it suppresses their density fluc- 
tuations and leads essentially to incompressibility of the 
vortex fluid. As demonstrated in the easy-plane quan- 
tum antiferromagnet on the triangular lattice ; 33 ' 34 in- 
cluding fluctuations about the flux-smeared mean field 
enables one to access a critical spin liquid with gap- 
less vortices. The theory has the structure of a (2+1)- 
dimensional [(2+l)D] quantum electrodynamics (QED3), 
with relativistic fermionic vortices minimally coupled to 
a non-compact U(l) gauge field. The resulting algebraic 
vortex liquid (AVL) phase is a novel critical spin liq- 
uid phase that exhibits neither magnetic nor any other 
symmetry-breaking order. This approach also allows one 
to study many competing orders in the vicinity of the 
gapless phase. 

When applied to the spin- 1/2 easy-plane antiferro- 
magnet on the kagome lattice, the duality transforma- 
tion combined with fcrmionization and flux smearing also 
leads to a low-energy effective QED3 theory with 8 flavors 
of Dirac fermions with an emergent SU(8) flavor sym- 
metry Amusingly, AVL's with SU(2)22, SU(4)22i24 and 
SU(6)^ 9 - emergent symmetries were obtained previously 
for quantum XY antiferromagnets on the triangular lat- 
tice, with integer spin, half-integer spin and half-integer 
spin in an applied magnetic field, respectively. QED 3 
theory with N flavors of Dirac fermions is known to re- 
alize a stable critical phase for sufficiently large N > N c . 
While numerical attempts to determine N c are so far 
inconclusive ; 40 ' 41 an estimate from the large-TV expan- 
sion suggests N c ~ 4i22. It seems very likely that N = 8 
is large enough, implying the presence of a stable criti- 
cal spin liquid ground state for the easy-plane spin-1/2 
quantum antiferromagnet on the kagome lattice. 

Although the algebraic vortex liquid and the algebraic 
spin liqui d 30 ' 31 obtained for the kagome lattice are ac- 
cessed in rather different ways, they share a number 
of commonalities, both theoretically and with regard to 
their experimental implications. Both approaches end 
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with a QED3 theory, the former a non-compact gauge 
field theory with fermionic vortices carrying an emergent 
SU(8) symmetry, and the latter a compact gauge the- 
ory with fermionic spinons with SU(4) symmetry. The 
non-compact nature of the gauge field in the AVL follows 
from the fact that the S z component of spin appears as 
the gauge flux in the dualized theory Hence the states 
with zero flux and 2tt flux are physically distinct; more- 
over, it follows that since the total S z is conserved, so also 
is the total gauge flux. This should be contrasted with 
the slave particle approach, where a compact gauge the- 
ory arises on the lattice and there is no such conservation 
law. Physically, this means that dynamical monopole op- 
erators, which could potentially destabilize the spin liq- 
uid and open a gap in the excitation spectrum, are not 
allowed in our low-energy theory for the AVL, while they 
are allowed in a low-energy description of spin liquids 
obtained using a slave particle framework. 

With regard to experimentally accessible quantities, 
both theories predict a power law specific heat C ~ T 2 
at low temperatures, which follows from the linear dis- 
persion of the fermions, and which would dominate the 
phonon contribution to the specific heat for magnets with 
exchange interactions significantly smaller that the De- 
bye frequency^ Because the spinons couple directly to 
magnetic fields, one would expect the specific heat in the 
ASL to be more sensitive to an external applied magnetic 
field than the AVL phase. Both theories predict a ther- 
mal conductivity which vanishes as k~ T" ) 46 ' 47 which is 
an interesting experimental signature reflecting the dy- 
namical mobility of the gaplcss excitations. 

The momentum resolved dynamical spin structure fac- 
tor which can be extracted from inelastic neutron ex- 
periments can in principle give very detailed information 
about the spin dynamics. Following the framework de- 
veloped on the triangular lattice ) 32 ' 33 i 34 ' 39 for the kagomc 
AVL we can extract the wavevectors in the magnetic Bril- 
louin zone which have gapless spin carrying excitations, 
and find 12 of them as shown in Fig.QJ In contrast, the 
kagome ASL phase is predicted to have gapless spin exci- 
tations at only 4 of these 12 momenta. The momentum 
space location of the gapless excitations is perhaps the 
best way to try and distinguish experimentally between 
these two (and any other) critical spin liquids. 

The rest of the paper is organized as follows. We start 
by introducing our model in Sec. [TTJ The low-energy ef- 
fective field theory is then developed in Sec. IIIII In Sec. 
IIV1 the properties of the critical spin liquid phase (AVL 
phase) are discussed, especially, the spin excitation (ro- 
ton spectrum) (Fig. [5]) , the S z dynamical structure fac- 
tor (Eq. ll9p . and the in-plane dynamical structure factor 
(Eq. [H]). We conclude in Sec. |Vl Some details of the 
analysis are presented in three Appendices. 



II. EXTENDED KAGOME MODEL 

Our primary interest is the kagome lattice spin-1/2 
model with easy-plane anisotropy (quantum XY model) 

n = ^[v^+^l+Ew. (!) 

r,r' r.r' 

with J r r , > J z r , . We consider the system with dominant 
nearest- neighbor exchange Ji, but will also allow some 
second- and third- neighbor exchanges J2 and J3. 

To set the stage, the model with only nearest-neighbor 
coupling has an extensive degeneracy of classical ground 
states, which is lifted by further- neighbor interactions. 
For example, with antiferromagnetic J2 > 0, the so- 
called g = phase is stabilized (Fig. [I}. On the other 
hand, for ferromagnetic J2 < 0, the so-called y/3 x \/3 
structure is the classical ground state (Fig.QJ. These or- 
dered phases are also realized in the quantum spin-1/2 
Heisenbcrg model for large enough J^, while the case 
with J2/J1 — is most challenging as mentioned in the 
introduction. 

The easy-plane spin system can be readily reformu- 
lated in terms of vortices. But in order to apply the 
fermionized vortex approach more simply, we consider a 
wider class of models which includes the nearest-neighbor 
kagomc model. Specifically, we add one extra site at the 
center of each hexagon, on which we put an integer spin 
(Fig. [T]). (Technical reasons for doing this are discussed 
at the end of Sec. IIIII ) In the rotor representation of 
spins, with a phase tp canonically conjugate to the inte- 
ger boson number, n ~ S z + 1/2 on the kagome sites and 
n ~ S z at the centers of each hexagon, such an extended 
model reads 

H = ^2 J r.r' cos((^ r - ip r >) + ^2 U r (n r - n°) 2 , (2) 

r.r' r 

where n° = 1/2 for the kagome lattice sites, while n° = 
for the added hexagon center sites. We have dropped the 
J z term for simplicity since it amounts to a mere renor- 
malization of the interactions in the effective field theory 
that we will derive. As shown in Fig.[TJ we take the inter- 
action to be J for the nearest neighbor half-integer spins, 
while the coupling that connects integer and half-integer 
spins is J'. When the on-site interaction U is infinitely 
large, U — > oo, there remain two low-energy states real- 
izing the Hilbert space of a S = 1/2 spin at the kagomc 
sites, whereas the frozen sector with n T = is selected at 
the hexagon center sites. In the quantum rotor model, 
we "soften" this constraint and take U to be finite. 

The additional integer-spin degrees of freedom do not 
spoil any symmetries of the original model, which for 
the record include lattice translations, reflections, and 
rotations, as well as XY spin symmetries. If we integrate 
over the extra degrees of freedom, the model looks very 
close to the original one. For example, in the limit when 
J' <C U, wc obtain the kagomc model with additional 
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FIG. 1: (Top) The quantum spin model with easy-plain 
anisotropy, Eq. ([2]), on the kagome lattice (thick line) sup- 
plemented by an extra site at the center of each hexagon. 
There is a half-integer spin for each kagome site whereas inte- 
ger spins are placed at each center of hexagons. Four spins in a 
unit cell are labeled by j = 1,2, 3, 4. For this triangular exten- 
sion of the kagome lattice, the dual lattice is the honeycomb 
lattice indicated in the upper left part. (Bottom) The q — 
state (Left) and the y/3x y/3 state (Right). Positive/negative 
vorticities (chiralities) for each triangle are denoted by +/— . 



ferromagnetic second and third neighbor exchanges J2 = 
J3 = —(J') 2 / (2^0, together with a renormalized nearest- 
neighbor exchange J\ = J — {J') 2 / (2f7). Such exchanges 
would slightly favor the -s/3 x V3 state in the classical 
limit, but as we will see, the spin liquid phase that we 
obtain by analyzing the extended kagome model will have 
enhanced correlations corresponding to both the q = 
and y/3 x y/3 phases that are nearby. Thus the small bias 
introduced by the added integer-spin sites appears to be 
not so important for the generic spin liquid phase that 
we want to describe. 



III. FERMIONIZED VORTEX DESCRIPTION 

The triangular extension of the kagome XY antiferro- 
magnet discussed above is somewhat similar to the tri- 
angular lattice spin-1/2 XY antiferromagnet studied in 
Refs. I33ll34l . except that one quarter of the sites is occu- 
pied by integer spins. The duality transformation for the 
present model proceeds identically to Refs. I33U341 . and 
turns the original spin Hamiltonian |(5J) into the dual 
Hamiltonian describing vortices hopping on the honey- 
comb lattice and interacting with a gauge field residing 



on the dual lattice links 



U 



(2 , )2 E( V >< a )r 
xx' r 

2 ^2 £ XX ' cos (0 X - 6» x / - a xx , - a xx ,) 



(3) 



Here x, x' label sites on the honeycomb lattice (see 



Fig.©, 



o ±i0 



is a vortex creation/annihilation operator 



at site x, and a xx ' , e xx ' represent a gauge field on the 
link connecting x and x'. J xx / is given by J XX ' = J(J') 
when the dual link (xx') crosses the original lattice link 
(rr') with J rr > = J(J'). The last term in the dual Hamil- 
tonian represents the vortex hopping where the hopping 
amplitude t xx / is given by f xx < = t(t') when the link (xx') 
crosses the link (rr') with J rr > = J (J')- Crudely, we have 
t/t' ~ J' j J since vortices hop more easily across weak 
links (Fig. [2]) . Finally, the static gauge field a encodes 
the average flux seen by the vortices when they move; 
this is described in more detail below. 

In the dual description, because of the frustration in 
the spin model, the average density of vortices is one half 
per site. On the other hand, the original boson density 
is viewed as a gauge flux 



E 



(4) 



(xx') 



Therefore, vortices experience on average 7r flux going 
around the triangular lattice sites with half-integer spins, 
but they see zero flux going around the sites with integer 
spins. Thus, one quarter of the hexagons will have zero 
flux as seen by the vortices, and this is where the present 
model departs from the considerations in Refs. [33.34. 

To treat the system of interacting vortices at finite den- 
sity, we focus on the two low-energy states with vortex 
number N x = and Y x = 1 at each site and view vortices 
as hard-core bosons. We can then employ the fermioniza- 
tion as in Refs. 133||34| and arrive at the following hopping 
Hamiltonian for fermionized vortices <i x 



Wf e 



kcx'dt^e-< a * 



,+a",+A x 



+ h.c. 



(5) 



where we have introduced a Chern-Simons field A whose 
flux is tied to the vortex density, (V x A) x = 2ttN x . 

Before proceeding with the analysis of the fermion- 
ized vortex Hamiltonian, we now point out the techni- 
cal reasons for considering the extended kagome model. 
If we apply the duality transformation to the kagome 
model with nearest-neighbor antiferromagnetic coupling 
Ji only, we would obtain a dual vortex theory on the 
dice lattice, which is the dual of the kagome lattice. In 
this theory, in order to reproduce the rich physics while 
restricting the vortex Hilbert space at each site to some- 
thing more manageable like the hard-core vortices de- 
scribed earlier, we would need to keep two low-energy 
states of vortices for each triangle of the kagome lattice 
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FIG. 2: The flux-smeared mean field background for fermion- 
ized vortices. Each hexagon is threaded by either zero or n- 
flux. (hexagons with zero-flux are specified by "0" whereas 
all the other hexagons are pierced by 7r-flux.) A convenient 
choice of a gauge is also shown: for links labeled by <-», we 
assign the Peierls phase factor e~ la *-*' — e I7r . The unit cell 
in this gauge consists of 16 sites as labeled in the figure. For 
links represented by a solid line we assign hopping amplitude 
t(= 1) whereas for links denoted by a broken line we assign 
hopping t' . We take x in the figure as an origin. 



(i.e., for each three-fold coordinated site of the dice lat- 
tice) , whereas there are three such states to keep for each 
hexagon (i.e., six- fold coordinated dice site). The lat- 
ter degrees of freedom are rather difficult to represent in 
terms of fermions. On the other hand, in the dual treat- 
ment of the extended model, such a six-fold coordinated 
dice lattice site is effectively split into six sites. Each 
such new site now has two low-energy states but the sites 
are coupled together, and this provides some caricature 
of the original important vortex states on the problem- 
atic dice sites. This representation now admits standard 
fermionization, and it was "found" in a rather natural 
way without any prior bias regarding how to treat the 
difficult dice sites with three important vortex states. 



A. Flux smearing mean field and low energy 
effective theory 

The flux attachment implemented above is formally 
exact, and there have been no approximations so far, 
although we have softened the discreteness constraints 
on the gauge fields. To arrive at a low-energy effective 
field theory, we now pick a saddle point configuration 
of the gauge field, incorporating a flux-smeared mean- 
field ansatz, and later will include fluctuations around 
it. That is, we first distribute (or "smear") the 2tt 
flux attached to each vortex on a dual site to the three 
(dual) plaquettes surrounding it. Since the vortices are at 
half filling and each hexagon contains six sites, the total 
smeared flux through each hexagon is 2tt, which is equiva- 
lent to no additional flux. Therefore we assume A xx ' = 
in the mean field. Also, we take a xx < = on average, in- 
cluding the fluctuations after we identify the important 
low-energy degrees of freedom for the fermionic vortices. 



Thus, at the mean field level, the original Hamiltonian 
^) is converted into a problem of fermions hopping on 
the honeycomb lattice, 



Hi, 



'erm.MF 



'did. 



h.c. 



(6) 



In the mean field, 3 out of 4 hexagons are pierced by ir- 
flux, while the remaining hexagons with an integer spin 
at the center have zero flux, see Fig. [5] Our choice of 
the gauge is shown in Fig. [2] The unit cell consists of 16 
sites, the area of which is twice as large as the physical 
unit cell of the original kagome spin system. 

This Hamiltonian can be solved easily by the Fourier 
transformation. The lattice vectors that translate the 
unit cell are given by 



Ei = (4,0), E 2 



1,V3). 



(7) 



and the first Brillouin zone (BZ) is specified as the 
Wigner-Seitz cell of the reciprocal lattice vectors 



( 5'27l)' 



G 2 = 



2tt 

"•VS 



(8) 



This is shown in Fig. [3J Note that, on the other hand, 
the physical unit cell is defined by EP hys = Ei/2 = (2, 0) , 
and EP hys = E 2 = (-1, \/3) • Thus the physical BZ is de- 
termined as the Wigner-Seitz cell of the reciprocal lattice 
vectors Gf ys = (tt, tt/V3) and Gf ys = (0, 2n/V3) , and 
is also shown in Fig. [3] and Fig. |4j 

Because of the particle-hole symmetry of the origi- 
nal quantum rotor Hamiltonian, which in terms of the 
fermionized vortices is realized as a vortex particle-hole 
symmetry (see Appendix [B]), the band structure of the 
flux-smeared mean-field Hamiltonian ([6]) is symmetric 
around zero energy, and consists of eight bands with 
positive energy and eight with negative energy. The 
7th and 8th bands have negative energy spectra and are 
completely degenerate^ each having four Fermi points 
Qi 2 3,4 that touch zero energy. In the same way, the 
9th and 10th bands have positive energy spectra and are 
completely degenerate, and each of them touches zero 
energy at the same four Fermi points Qi.2,3,4 thus com- 
pleting the Dirac nodal spectrum. With the gauge choice 
specified in Fig. [2J the locations of the Fermi points in 
the BZ are 



, 7T 7T \ 



Q 2 = Q 4 



/ 7T 377 \ 

T2'~47!y ' 



(9) 
(10) 



(see Fig. [3]). The low-energy spectrum consists of eight 
gapless two-component Dirac fermions with identical 
Fermi velocities given by (here we take t = 1 for con- 
venience) 



v F 



(3 + i' 2 )(3 + 2*' 2 ) 



(11) 
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FIG. 3: The first Brillouin zone for the flux smeared mean 
field ansatz and the location of the nodes Qi,...,4. The hexag- 
onal physical Brillouin zone is also presented. Small hexagons 
are guide for eyes. 



Hence the low-energy effective theory enjoys an emer- 
gent SU(8) symmetry among eight Dirac cones, which is 
protected at the kinetic-energy level by the underlying 
discrete symmetries in the problem (see Appendix [Cjl . 

The lattice fermionized vortex operator is conveniently 
expanded in terms of slowly varying continuum fields 

^iaa( x ) b y 



EE 

ct=l i=I,II o=l,2 



E 



lka«(*). (12) 



where indices a = 1, . . . , 4, % = I, II, and a = 1,2, refer 
to the four Fermi momenta, two sublattices of the hon- 
eycomb lattice, and two degenerate bands, respectively, 
while n = 1, . . . , 16 specifics the site label in the unit 
cell, cf. Fig. [2] The 16-component wavefunctions 4> iaa (n) 
represent modes at the Fermi points and are specified in 
Appendix [Al 

Working with these continuum fields, we find that the 
low-energy Hamiltonian has the same form vf(— <J X p x + 
u 2 p y ) near each Dirac node, where a 12,3 are Pauli ma- 
trices that act on the "Lorentz" indices i = I, II. Cor- 
respondingly, by introducing the gamma matrices in 
(2 + l)D 

7o = cr 3 , 71 = -cr 2 , 72 = -er 1 , 
7/^+7^ = 2^, fi,v = 0,1,2, (13) 

the Euclidean low-energy effective action is given by 



f 8 

S eS = J d 3 x 22 tp (a) 7/AVv 



(14) 



(«)=! 



where (a) — aa is a collective index running from 1 to 
8; x^ =0 ' 1,2 represents imaginary time as well as spatial 
coordinates; — "0^7o; and we have scaled out vp. 



Having identified the low-energy fcrmionic degrees of 
freedom, we now reinstate gauge field fluctuations. The 
complete effective action is that of SU(8) QED3 coupled 
in addition with the Chern-Simons field A and given by 



S e f[ = I d x 



(«)=! 
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Here we also included possible fcrmion bilincars C2 / and 
four fermion interactions £4/. e 2 ~ U" 1 stands for the 
coupling constant of the gaug e field a^. 

As discussed in Refs.[HIU the Chern-Simons field A^ 
is irrelevant. The microscopic symmetries of the spin 
Hamiltonian prohibit fcrmion mass terms £2/ from ap- 
pearing in the effective action (sec Appcndix[C|). Further- 
more, four fermion interactions C^f in Eq. (|15p arc irrele- 
vant when the number of fcrmion flavors is large enough, 
N > N c . Thus, provided N c < 8, we have obtained a 
stable Algebraic Vortex Liquid phase that is described 
by QED3 with an emergent SU(8) symmetry. Using the 
available theoretical understanding of such theories, we 
will now describe the main properties of the AVL phase, 
which in terms of the original spins is a gapless spin liquid 
phase. 



IV. PROPERTIES OF THE AVL PHASE 
A. Roton spectrum 

First note that the specific momenta of the Dirac 
points Eq. (JTUJ) and Fig. [3] are gauge-dependent. We are 
interested in the physical properties of the system, which 
are gauge invariant. One such property is the vortex- 
antivortex excitation spectrum: By moving a vortex from 
an occupied state k with energy E~ (k) to an unoccupied 
state k' with energy £ ,+ (k'), we are creating an excita- 
tion that carries energy E + (k') — E~(\c). Some care is 
needed if we want to construct such a state with a definite 
momentum in the physical BZ, since we need to consider 
both k' — k and k' k + Gi. By examining Fig. [3J we 
find that the vortex-antivortex continuum goes to zero 
at twelve points 0, ±Q, Mi j2 ,3, ±Pi,2,3 in the BZ of the 
kagome lattice shown in Fig. [U 

An accompanying plot of the lower edge of the vortex- 
antivortex continuum, 

A rot (p) = min k {i;+(k + p)-£-(k)}, (16) 

which can be interpreted as a roton excitation energy, is 
shown in Fig. [5] (top panel) along several BZ cuts. Var- 
ious proximate phases, e.g., magnetically ordered states, 
will instead have gapped vortices. But if such a gap is 



small, we expect that the lower edge of the full excita- 
tion spectrum will be dominated by deep roton minima 
near the same wavevectors (except where the rotons are 
masked by spin waves). 

As a contrasting example, in Fig. [5] (bottom panel) we 
also show the spin-wave spectrum in the q = magneti- 
cally ordered phas o 49 ' 50 , assuming this state is stabilized 
by some means. (Actually, starting from any of the clas- 
sically degenerate ground states of the nearest-neighbor 
kagome XY model leads to the same linear spin-wave 
theory if one rotates the spin quantization axis suitably 
for each site.) There are three branches. For the XY 
antifcrromagnct, these are 

w$ )3 (p) = 5J v /4-A 1 , 2 , 3 (p), 
Ai(p) = -2, 
A 2 , 3 (p) = IT \/l + 8cospi cosp 2 cosp 3 , (17) 

where p x = p x , p 2 = p x /2 + y/3p y /2, and p 3 = p x /2 - 
s/3p y /2. Clearly, the spin-wave frequency is large and of 
order J throughout the BZ except near zero momentum. 
It is possible that higher order in 1/5 spin wave correc- 
tions will lead to roton minima. It would be interesting 
to explore the possibility of such roton minima in the 
magnetically ordered phase to see if they are coincident 
with the gapless momenta of the AVL phase as shown 
in Fig. 3J Such a coincidence between the roton minima 
obtained with spin wave perturbation theor y 36 i 37 ' 38 and 
with the AVL approach 3 ^ was indeed found within the 
magnetically ordered phase of the square lattice spin- 1/2 
quantum antifcrromagnet with a weak second neighbor 
exchange across one diagonal of each elementary square 
(a spatially anisotropic triangular lattice). It would be 
useful to look for such signatures in the kagome experi- 
ments. 

The momenta of the low-energy or gapless excita- 
tions capture important lattice- to intermediate-scale 
physics and provide one of the simplest characteristics 
that can be used to distinguish between different spin 
liquid phases. For example, in the algebraic spin liquid 
state proposed by Hastings and Ran et aZi 30 i 31 for the 
5 = 1/2 nearest-neighbor Heisenberg model, the gapless 
spinon-antispinon excitations occur at four momenta {0, 
Mi ; 2,3} in the BZ. On the other hand, Sachde\*2i and 
Wang and Vishwanath^ discuss four different gapped Z 2 
spin liquid phases on the kagome lattice using Schwingcr 
bosons. One of the states has minimum in the gapped 
spin excitations at 0; one has minima at ±Q; and the 
two remaining spin liquids have minima at the same 12 
momenta as the AVL phase, Fig. [4] 

B. S Z S Z correlator 

A more formal approach towards characterizing the 
system with gauge interactions is to deduce the PSG 
transformations for a particular mean-field stated 
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y 




FIG. 4: Summary of the main characterizations of the 
kagome AVL phase. The low-energy excitations are lo- 
cated at wavevectors 0, ±Q, Mi : 2,3, and ±Pi,2,3 in the 
physical Brillouin zone of the kagome lattice, where Q = 
(*r/3,7r/>/3); Mi = (77/2,-77/(2^)), M 2 = (77/2,77/(2^)), 
M 3 = (0,77/^3); Pi = (tt/3,0), P 2 = (-77/6,77/(2^3)), and 
P3 = (-77/6, -77/(2^)). Both S z and S + exhibit power law 
correlations at all these momenta, Eqs. (|19|l and (|22[) . The S z 
correlations are "enhanced" {n z ~ 2.46) for the subset ±Q, 
±Pi,2,3 and not enhanced (rj z = 3) for the remaining wavevec- 
tors. On the other hand, all S + correlations are characterized 
by the same exponent rj± ~ 3.24. 




r m k r 

(Px» Py) 

FIG. 5: (Top) Lower edge of the vortex-antivortex continuum 
along several cuts in the Brillouin zone as shown in Fig. [4] for 
parameter values t'/t = 0.5 and t'/t = 0.2 from the higher 
to lower curve. (Bottom) Spin wave excitations, Eq. (|17[) . of 
the kagome XY antiferromagnet with the nearest-neighbour 
exchange J only along the same cuts in the Brillouin zone. 
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Specifically, the notion of symmetry is enlarged to also 
include gauge transformations in order to maintain in- 
variancc of the mean-field Hamiltonian under the usual 
symmetry operations. 

Some details are presented in Appendix [Cj where we 
show the calculated transformation properties of the con- 
tinuum fcrmion fields ?/>, ip. In particular, with the PSG 
analysis, the low-energy roton excitations can be studied 
by analyzing fermion bilinears. 

The S^-componcnt of the spin is expressed as a gauge 
flux in our QED3 theory of vortices. Both the flux in- 
duced by the dynamical U(l) gauge field itself and the 
flux induced by vortex currents contribute to the S Z S Z 
correlation. Formally, the S z operator can be expressed 
in terms of the gauge flux and appropriate fcrmion bilin- 
ears ipQip representing the contributions by vortex cur- 
rents, 

Sr ~ ^ + E A ^ Kl ' r ^+-" ( 18 ) 



where the summation over i runs over bilinears ipGiip 
with appropriate symmetry properties so that the right 
hand side transforms identically to S z . The factor e ~ %v " vV 
represents the momenta carried by these bilinears, which 
are the same as the roton minima shown in Fig.|H and Ai 
is some amplitude for the bilinear Qi [in fact, one needs 
to specify four such amplitudes Ai(J) since there are four 
sites j = 1, ... ,4 in the unit cell, Fig. [T]. 

Assuming the N = 8 QED3 is critical, all correlation 
functions decay algebraically. We then expect to find 
power law correlation at all momenta shown in Fig.|U As 
a consequence, the singular part of the structure factor 
at zero temperature near such wavevector K is given by 

Q(uj 2 - a 2 ) 

S zz (k = K + q, U ) cx ( ^ q2)(2 ^ )/2 . (19) 

The exponent r\ z characterizes the strength of the S z 
correlations, and smaller values correspond to more sin- 
gular and therefore more pronounced correlations. In 
QED 3 , the corresponding exponent value for the flux 
V x a that contributes to the S z correlation near zero 
momentum is r/ z = 3, and the same exponent also char- 
acterizes fermionic bilinears that are conserved currents, 
see Appendix [O On the other hand, bilinears that are 
not conserved currents have their scaling dimensions en- 
hanced by the gauge field fluctuations as compared with 
the mean field. Upon analyzing the transformation prop- 
erties of the fermionic bilinears, we conclude that there 
are indeed such enhanced contributions to S z near the 
momenta ±Q and ±Pi :2 ,3 (explicit expressions are given 
in Appendix |C|) . The new exponent can be estimated 
from the large-TV treatments of the QED ^ 43 ! 44 to be 
r] z w 3 - 128/(37r 2 iV) = 2.46 for N = 8. 



C. S + S~ correlator 

Since S z ~ (V x a)/(2n), the spin raising operator S + 
in the dual theory is realized as an operator that creates 
2tt gauge flux. Following Ref.[13, we treat such flux inser- 
tion ( "monopolc insertion" ) classically as a change of the 
background gauge field configuration felt by fermions. By 
determining quantum numbers carried by the monopolc 
insertion operators, we can find the wavevectors of the 
dominant S + correlations. 

Here we only outline the procedure (for details, see 
Refs. I32ll34h . In the presence of 2tt background flux 
inserted at x = 0, each Dirac fermion species (a) = 
1 ... 8 has a quasi-localized zero-energy state near the in- 
serted flux. We will denote the creation operators for 
such fermionic zero modes by f) a y The correspond- 
ing two-component wavefunctions have the form ip(x) ~ 

T*T (, ) ' w ^ cn auows us to relate the transformation 

properties of /J , to those of the first component of 
(the latter are summarized in Appendix [B| . 

The monopole creation operators A4' are defined by 
the combination of the flux insertion and the subsequent 
filling of four out of eight fermionic zero modes in order 
to ensure gauge neutrality: 

where \DSq) and \DS+) are the Fermi-Dirac sea in the 
absence and presence of +2ir flux, respectively. There 
are total of 70 distinct such monopolcs. The quan- 
tum numbers of the monopole operators are then deter- 
mined by the transformation properties of the operators 

/(V(V(Vc<5) and those of \ DS +) relative t0 \ DS o)- 

Both \DS + ) and \DSo) arc expected to be eigenstatcs 
of the translation operator Tg r , and we denote the ra- 
tio of the eigenvalues as e tn ° ffs,rt ' 5r . On the other hand, 
we can diagonalize the action of T St on / ( t Q) / ( t /3) / ( t y) /( t 5) 

to get a set of eigenvalues e lUi ' 5r , i = 1, . . . , 70. Thus, 
the momenta of the monopole operators (at zero energy) 
arc determined by {11^} which are offset by the hitherto 
unknown wavevector n o ff sc t. However, once we exam- 
ine {IIj;}, we find that the offset is in fact uniquely fixed 
by the requirement that the monopole momenta are dis- 
tributed in the physical BZ in a way that respects all 
lattice symmetries. 

Monopolc momenta determined in this manner are 
found to be given by the very same wavevectors 0, ±Q, 
Mi,2,3, and ±Pi j2; 3, as the leading S z correlation (Fig. 
[4]). For the record, we find the following monopole mul- 
tiplicities given in the parentheses: (12 monopoles), 
±Q (5 each), Mi >2 ,3 (8 each), and ±Pi,2,3 (4 each). 

The monopole insertions which add S z — 1 have over- 
lap with the exact lowest energy spin carrying eigenstates 
in the AVL. At the lowest energies such eigenstates will 
be centered around the monopolc wavevectors in the Bril- 
louin zone, and we expect that they will disperse with an 
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energy growing linearly with small deviations in the mo- 
menta. But presumably any local operator that adds 
spin one will create a linearly dispersing but ovcrdamped 
"particle" excitation. More formally, we can expand the 
S + operator in terms of the continuum fields as follows, 



(21) 



where 11,; is the momentum of the monopole Mi. The 
right hand side has correct transformation properties un- 
der the translations, but not all of the monopole opera- 
tors contribute to S + when other symmetries are taken 
into account. For example, some of the zero momentum 
monopoles have opposite eigenvalues under the lattice in- 
version, and the same happens for the momenta Mi. 2,3- 
However, given the large monopole multiplicities, it ap- 
pears very likely that all momenta shown in Fig. [4] will 
be present in S + (but we have not performed an analysis 
of all quantum numbers so far). 

The in-plane dynamical structure factor S~* around 
each Hi is thus expected to have the same critical form 
as S zz , 



S + -{\t = n + q,w) oc 



0(u 2 - q 2 ) 
(cu 2 - q 2 )(2-')±)/2 ' 



(22) 



but the characteristic exponent rj± is now that appro- 
priate for the monopole operators and is expected to be 
the same for all momenta. Ref. Hi] calculated this in the 
large-TV QED 3 : r] mon w 0.53iV - 1. Setting N = 8, we 
estimate i]± f» 3.24, which is larger than even the non- 
enhanced exponent r] z = 3 entering the S z correlations. 

The above results with r]± > rj z suggest that within 
the AVL phase the system is "closer" to an Ising order- 
ing of S z than of an in-plane XY ordering of S + . This 
is rather puzzling for an easy-plane spin model, where 
one would expect XY order to develop more readily than 
Ising order. This conclusion is reminiscent of the classi- 
cal kagome spin model, where the easy-plane anisotropy 
present in the XY model destroys the zero temperature 
coplanar order of the Heisenberg model. However, this 
interpretation of the inequality r)± > r) z is perhaps some- 
what misleading. Indeed, as we discuss in Appendix [Cl 
condensation of any of the enhanced fermionic bilinears 
that contribute to S z would drive XY order in addition 
to S z order - forming a supersolid phase of the bosonic 
spins. Thus, care must be taken when translating the 
long-distance behavior of the AVL correlators into order- 
ing tendencies. 

The exact diagonalization study by Sindzingrc 35 finds 
that in the nearest-neighbor kagome XY antifcrromag- 
net the S + correlations arc larger than the S z correla- 
tions, which is different from the AVL prediction. It 
would be interesting to examine this further, perhaps 
also in a model with further-neighbor exchanges. With 
regards to such numerical studies, we also want to point 
out that the discussed dynamical spin structure factor, 
Eqs. (|19p and l[2"2"|). translate into equal-time spin corre- 
lations that decay as r -1-1 ' 2 -*. For the estimated values 



°f Vz,±i the decay is rather quick, and could be hard to 
distinguish from short-range correlations. 



V. CONCLUSION 

Using fermionized vortices, we have studied the 
easy-plane spin-1/2 Heisenberg antifcrromagnet on the 
kagome lattice and accessed a novel gapless spin liquid 
phase. The effective field theory for the resulting Al- 
gebraic Vortex Liquid phase, is (2+l)D QED 3 with an 
emergent SU(8) flavor symmetry. This large number of 
flavors can be thought to have its origin in many compet- 
ing ordered states that are intrinsic to the "loose" net- 
work of the corner-sharing triangles in the kagome lattice. 
It is also likely that for this number of Dirac fermions the 
critical QED3 theory is stable against the dynamical gap 
formation, which allows us to expect that the AVL phase 
is stable. 

The gapless nature of the AVL phase has important 
thermodynamic consequences. For example, we predict 
the specific heat to behave as C ~ T 2 at low tempera- 
tures. Since vortices carry no spin, we expect this con- 
tribution to be unaffected by the application of a mag- 
netic field. Interestingly, such behavior was observed in 
SrCr8_; C Ga4+ x Oi9 and was interpreted in terms of sin- 
glets dominating the many-body spectrunii&ii Since the 
vortices are mobile and carry energy, we furthermore pre- 
dict significant thermal conductivity, which would be in- 
teresting to measure in such candidate gapless spin liq- 
uids. 

The more detailed properties of the kagome AVL phase 
were studied by the PSG analysis, leading to detailed pre- 
dictions for the spin structure factor. We found that the 
dominant low energy spin correlations occur at twelve 
specific locations in the BZ shown in Fig. 3J These 
wavevectors encode important intermediate-scale physics 
and can be looked for in experiments and numerical stud- 
ies. They may be also used to compare and contrast with 
other theoretical proposals of spin liquid states. 

The recent experiment a 19 i 20 ' 21 for the spin-1/2 kagome 
material ZnCu3(OH)gCl2 observed power-law behavior in 
lo of the structure factor S(k, uf) in a powder sample. The 
specific heat measurements and local magnetic probes 
also suggest gapless spin liquid physics in this compound. 
Unfortunately, no detailed momentum-resolved informa- 
tion is available so far. It should be also noted that 
our results were derived assuming easy-plane character 
of the spin interactions, and do not apply directly if 
the material has no significant such spin anisotropy^ 
Hopefully, more experiments in the near future will pro- 
vide detailed microscopic characterization of this mate- 
rial, such as the presence/absence of the Dzyaloshinskii- 
Moriya interaction^, and clarify the appropriate spin 
modeling. In the case of the Heisenber g spin symme- 
try and nearest-neighbor exchanges, Refs. I3dll3ll offer an- 
other candidate critical spin liquid for the kagome lattice, 
obtained using a slave fermion construction. This algc- 
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braic spin liquid (ASL) state differs significantly from the 
presented AVL phase as discussed in some detail at the 
end of the introduction and towards the end of Sec. lIV Al 
Briefly, within the ASL the momentum space locations of 
the gapless spin carrying excitations are only a small sub- 
set of those predicted for the AVL. Moreover, while both 
phases support gapless Dirac fcrmions that contribute a 
T 2 specific heat, in the AVL phase the fermions (vor- 
tices) are spinless in contrast to the fcrmionic spinons in 
the ASL. Thus, one would expect the specific heat to be 
more sensitive to an applied magnetic field in the ASL 
than in the AVL. 

Contrary to SrCr 8 _ 2 .Ga4 +a .Oig, the experimental data 
on ZnCu3(OH)eCl2 shows that the specific heat is af- 
fected by the magnetic field in the low temperature 
regime^. However, as pointed out by Ran et al.^ 
T ~ 10K is likely to be an appropriate temperature scale 
to test several theoretical approaches to this material, 
since spin liquid physics might be masked by, for exam- 
ple, impurity effects, Dzyaloshinskii-Moriya coupling, or 
other complicating spin interactions, at the lowest tem- 
perature scales. (Indeed, the susceptibility data are per- 
haps consistent with the presence of impurities, although 
it is unclear at this stage if the peculiar temperature de- 
pendence of the specific heat also has its origins in impu- 
rities.) The theoretical prediction of the T 2 specific heat 
is consistent with the experiments for T > 10K. 

A general outstanding issue is the connection, if any, 
between the critical spin liquids obtained with slave 
fermions and with fermionic vortices. Perhaps our AVL 
phase corresponds to some algebraic spin liquid ansatz, 
but at the moment this is unclear. 
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APPENDIX A: ZERO MODE WAVE FUNCTIONS 
AT THE FERMI POINTS 

In this appendix, the 16-component wavefunctions 
4>i a a( n ) (n = 1, . . . , 16) representing zero modes at the 
Fermi points are presented explicitly. These wave func- 
tions are necessary for the PSG analysis. 

At the node Qi, 



at the node Q2 



'11,1,1 



.k/a I Xa+ 




Qi 

in/12 ( Xa- 
^ 



,— iw/4 



Cjl I ) VII, 1,2 
Xb + 



a -wr/12 



Qi 



(Al) 
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at the node Q 3 , 
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and at the node Q4, 



"1,4,1 



"11,4,1 





Q4 

Xb+ 



J 1A,2 



'11,4,2 = e 



1 xt 

( X*A- \ 



(A4) 



where a convenient choice for the eight-component zero 
energy wave functions xa and xb is 



xts = Na 



1 

a 



\ 



Xts = N B 



(xy) 2 (y 4 + (xy*) 2 ) 
-a{xy) 2 (y 4 - {xy*) 2 ) 
-xy* 3 + axy*(y A - {xy*) 2 ) 
H [x*y* 3 + axy* 3 ] 
-~x 3 y* (y 4 + (xy*) 2 ) - ax*y* 3 
V t'x 3 y* [y 4 x* 2 +y* 2 -a (y 4 - {xy*) 2 )] J 

(A5) 



/ 1 \ 

P 

(x*y*) 2 {y* 4 + (x*y) 2 ) 
-f3(x*y*) 2 {y* 4 -(x*y) 2 ) 
-f [x*y 3 + f3xy 3 ] 
x*y(y* 4 + {x*y) 2 ) + /3x*y 3 
t'x* 3 y [(3 (y* 4 x 2 - y 2 ) - y* 4 - {x*y) 

xy 3 - [3x* 3 y (y* 4 - (x*y) 2 ) J 



Here x := e lfc */ 2 , y := e lk y/^ 2 ^ for Q = (k x , k y ), 

s3V2x* 3 y* 3 - (1 + x e y 6 ){x 4 + x* 4 ) 
a ~ 2x- 2 (3 - x 2 y* e - x* 2 y 6 ) ' 

si3V2x* 3 y* 3 - (1 - x 6 y 6 ){x 4 + x* 4 ) 



(A6) 







2x- 2 {3 - x 2 y* 6 - x* 2 y 6 ) 



, (A7) 



and Na and Ng are t'-dependent normalization factors. 
These zero modes are chosen in such a way that the s = 
±1 wavefunctions arc orthogonal to each other for any t' . 
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In order to work with the continuum Dirac fields 
Eq. (O, it is convenient to introduce four sets of the 
Pauli matrices which act on different gradings. In the 
following, Pauli matrices denoted by p K (re = 0,1,2,3) 
act on Latin indices a = 1,2 (p° is an identity ma- 
trix). Pauli matrices denoted by a K (and hence 7 K in- 
troduced in Eq. [T3")) act on Latin indices i = I, II. We 
separate the node momenta into two groups (Qi,Qa) 
and (Q2,Q4)- Then, Pauli matrices denoted by r K act 
within a given group, either (Qi,Qa) or (Q2,Q4)- Fi- 
nally, Pauli matrices denoted by p K act on the space 
spanned by (Qi3,Q2 4)- We will also use the notation 
r± = r 1 ± ir 2 . 



APPENDIX B: PROJECTIVE SYMMETRY 
GROUP (PSG) ANALYSIS 

Any symmetry of the original lattice spin model has a 
representation in terms of vortices. Unfortunately, upon 
fermionization, time-reversal symmetry transformation 
becomes highly non-local and we do not know how to 
implement it in the low-energy effective field theory. 

Since a specific configuration of the Chcrn-Simons 
gauge field was picked in the flux smearing mean field, we 
had to fix the gauge and work with the enlarged unit cell. 
The spatial symmetries of the original lattice spin model 
are, however, still maintained if the symmetry transfor- 
mations are followed by subsequent gauge transforma- 
tions - the original symmetries become "projective sym- 
metries" in the effective field theory. The PSG analysis 
is necessary when we try to make a connection between 
the original spin model and the effective field theory. 

The original spin system has the following symmetries 
besides the global U(l): 

^Ei/2 : translation by Ei/2, 
Tg 2 : translation by E2, 

: rotation by 7r around a kagome site, 
1Z X : reflection with respect to y axis, 
C : particle-hole, 

T : time-reversal. (Bl) 

We can also include a rotation by 7r/3 about honeycomb 
center but will not consider it here. 

The first four transformations Tg r , Re-, 1Z X act on spa- 
tial coordinates r = (x, y), r — > r + <5r, r — > (xcosO — 
ysm0,xsm0 + ycos8), r — > (— x,y), respectively. On the 
other hand, the particle-hole symmetry C and the time 
reversal T act on spin operators as 

C: S z -► ~S Z , S x ± iS y -> S x T iS y , 

T : S^-S, (B2) 

Here, the anti-unitary nature of the time reversal is re- 
flected in its action on the complex number i — > — i. 



Symmetries in terms of the rotor representation can 
then be deduced as 

C : n — > —n, tp — ► —95, 

T : n — > — n, <p — > tp + 7r, i — > —i. (B3) 

Symmetry properties of bosonic and fermionic vor- 
tices are deduced from their defining relations. Due 
to the Chern-Simons flux attachment, the mirror and 
time-reversal symmetries are implemented in a non-local 
fashion in the fermionized theory. However, if we com- 
bine T with 1Z X and C, the resulting modified reflection 
7Z X := 1Z X CT can still be realized locally. We summarize 
the symmetry properties of the fermions d, dJ in Table [U 
We also introduce a formal fermion time reversal by 

Tfcrm : d ^ d, i — > -i. (B4) 

The necessary explanations are the same as in Refs.llH 
and are not repeated here. 

Finally, the symmetry properties of the slowly varying 
continuum fields can be deduced from Eq. (TT2")) . As 
explained, these are realized projectively^ and the sym- 
metry transformations for the continuum fermion fields 
are summarized in Table [TT1 

APPENDIX C: FERMION BILINEARS 

Once we determine the symmetry properties of the 
continuum fermion fields i/},ip, we can discuss symmetry 
properties of the gauge invariant bilinears tpQtp where Q is 
a 16 x 16 matrix. Using the gradings specified at the end 
of Appendix El the bilinears can be conveniently written 
as 

?K> a t q /V ( C1 ) 

where re, A, a, (3 run from to 3 and hence there are 
4 4 = 64 x 4 such bilinears. Among these bilinears, 64 x 3 
bilinears with the Lorentz index re = 1, 2, 3 represent U(8) 
currents, 

V^^V^VV (C2) 

On the other hand, there are 64 remaining bilinears with 
re = 0, 

$a°p X T a pPip. (C3) 

Since the bilinears (|C2[) are conserved currents, they 
maintain their engineering scaling dimensions even in the 
presence of the gauge fluctuations. On the other hand, 
the bilinears of the form (|C3[) ( "enhanced bilinears" ) de- 
velop anomalous dimensions when we include the gauge 
fluctuations. 

Let us give some relevant examples of bilinears. The 
•\/3x vo magnetically ordered state (shown in Fig.[T]) cor- 
responds to the staggered charge density wave of vortices, 
which is obtained by adding the corresponding staggered 
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spin S 


rotor <p, n 


c 


S y _ S y ; S z _ S Z 


n — > — n, v 3 — * — V 3 


T 


S — > — S, i — > —i 


n — » —71, <y3 — > 1/3 + 7T, i — > — i 




S * S , 7 > 7, X ► — X 


7i — > +71, 95 — > — </? + 7r, i — > —i, a; — > —x 





vortex 8,N and gauge field a, e 


fermion d, d) 


c 


a -> -a, 6 -> -0, e -> -e, JV -» 1 - N 


d^ (-l)*d T 


T 


a — > —a, # — > — 0, e ^ e, N ^ N , i ^ —i 




ft* 


a — » —a, — » — 0, e — *• e, JV — > JV, i — > — i, x — > — x 


d — > d, a; — > — x, 7^—7 



TABLE I: Summary of symmetry transformations for spin, rotor, vortex, and fermionized vortex operators on the lattice. The 
sign factor (— l) 1 in the action of the particle-hole transformation on the fermionized vortices, C : d — > (— l) l d^, is 1 on one of 
the sublattice of the dual lattice whereas it is —1 on the other. 





T El/2 


Te 2 


ft* 




C 




4 - 






£7 3 p 2 T l e ^^/12 /i 3 v , 


tV [V^]"' 


2 2 / 
T (T If) 



TABLE II: Summary of symmetry transformations for continuum fermion fields (upto unimportant U(l) phase factors). 



chemical potential (i.e., ±Sp on the up/down kagome tri- 
angles). Using the listed symmetries, it is readily verified 
that the following two bilincars in the expression 

B V3 x Vs = a ^° ° A V V- + P i>a 3 fi 3 T°p 3 ^ (C4) 

transform in the manner expected of the staggered chem- 
ical potential (the coefficients a and (3 are generically in- 
dependent). This can be also checked explicitly by writ- 
ing such chemical potential in terms of the continuum 
fields obtaining some t'-dependent coefficients a and (3. 
In the similar vein, the q = state (also shown in Fig. [T]) 
corresponds to the uniform chemical potential for vortices 
on the up and down kagome triangles, which is realized 
with the following bilinears 

B q=0 = a ^a /j, 3 T 3 p 3 ip + V5cr 3 /r°pV • (C5) 

Note that each B^ x ^ and B q= o contains an enhanced 
bilinear of the QED3 theory, and therefore both orders 
are "present" in the AVL phase as enhanced critical fluc- 
tuations. It is in such sense that the enhanced bilinears 
encode the potential nearby orders. 

With the table[H]at hand, it is a simple matter to check 
that spatial symmetries (translation, rotations, and re- 
flection) and particle-hole symmetry prohibit all fermion 
bilinears from appearing in the effective action (|15p . ex- 
cept ipifj an d V'cr 3 p 3 T 3 p 3 ip. If we are allowed to require 
the invariance under 7f crm , then these bilincars would be 
prohibited as well. However, even if we do not use 7f erm , 
we exclude these bilinears from the continuum theory us- 
ing the following argument from Ref. l34l : Consider first 
adding the bilinear ipi/j to the action and analyze the re- 
sulting phase for the original spin model. This bilinear 
opens a gap in the fermion spectrum, and proceeding as 
in Ref. |3J we conclude that this phase is in fact a chiral 
spin liquid that breaks the physical time reversal. There- 
fore, if we are interested in a time-reversal invariant spin 



liquid, we are to exclude this term. The situation with 
the ^a 3 p 3 T 3 p 3 4> term is less clear since by itself it would 
lead to small Fermi pockets, and it is then difficult to 
deduce the physical state of the original spin system. In 
the presence of both terms, depending on their relative 
magnitude one may have either a gap or Fermi pockets. 
However, we can plausibly argue that these pockets tend 
to be unstable towards a gapped phase that is continu- 
ously connected to the same chiral phase obtained when 
the mass term tpip dominates. Since we are primarily in- 
terested in the states that are not chiral spin liquid (e.g., 
states that appear from the AVL description by sponta- 
neously generating some other mass terms like B^ x ^ 
or -B g= o)j wc drop the bilinears ipip an( i i><J 3 p?r 3 p 3 ^) 
from further considerations. The validity of this assump- 
tion and the closely related issue of neglecting irrelevant 
higher-derivative Chern-Simons terms from the final AVL 
action are the main unresolved questions about the AVL 
approach (see Ref. HH for some discussion). 

As another example of the application of the derived 
PSG, we write explicitly combinations of enhanced bilin- 
ears that contribute to S| ~ (Vxa)j/ (2n)+Gj+Fj + . . . : 

G.j = g 3 (e lQ ^e- l1T/l2 B+ + h.c.) , (C6) 
F,j = // fe^e 8 */ 1 '^ +h.c.) 

+fUe iP '-''e*"' 1 *B+ i +h.c.) 

+ff (e irVr V 57r/12 B+ +h.c.) . (C7) 

Here j refers to the "basis" labels in the unit cell con- 
sisting of four sites in the extended model shown in 
Fig. [TJ the wave vectors are the ones shown in Fig. [4] 
while the corresponding bilinears are Bq = ipp 2 T + p°tp, 

S+ = Tjjp 3 T+p 2 ij, B+ 2 = ijp 2 T+p 3 iP, B+ 3 = 4>p°T + p 2 ip. 
The above is to be interpreted as an expansion of the mi- 
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croscopic operators S* defined on the original spin lattice 
sites in terms of the continuum fields in the theory. It is 
straightforward but tedious to verify that one can choose 
real parameters <?j and /j 1 ' 2 ' 3 so that G and F have iden- 



tical transformation properties with V x a and therefore 
indeed contribute to S z . We thus conclude that S z has 
enhanced correlations at the wavevectors ±Q, ±Pi,2,3 as 
claimed in the main body. 
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